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Abstract 

We introduce an interesting hierarchy of rational order chaotic maps that posses an 
invariant measure. In contrast to the previously introduced hierarchy of chaotic maps 
[1, 2, 3, 4, 5], with merely entropy production, the rational order chaotic maps can simul- 
taneously produce and consume entropy . We compute the Kolmogorov-Sinai entropy 
of theses maps analytically and also their Lyapunov exponent numerically, where that 
obtained numerical results support the analytical calculations. 

Keywords: entropy production and entropy consumption, chaotic maps, chaos, 
Lyapunov exponent, Kolmogorov-Sinai entropy. 
PACs numbers:05.45.Ra, 05.45.Jn, 05.45.Tp 
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There has been some attempts [1, 2, 3, 4, 5, 6] at introducing the hierarchy of chaotic maps 
with an invariant measure in the recent years. The objective of these papers is to describe 
the dynamic behavior of chaotic maps using Kolmogorov-Sinai entropy. These hierarchies of 
chaotic maps are of interest as models for describing of behavior of dynamical systems. As 
an example, the random chaotic maps have attracted the attention of physicists as models 
of convection by temporarily irregular fluid flows [7]. Once a map is determined, the long 
term statistical behavior is described by a probability density function, which can be obtained 
either by solving the Frobenius-Perron equation [8] or can be estimated by measurement of the 
system. Therefore, the complexity, non-linearity and non-stationarity of physical, chemical, 
biological, physiological and financial systems [9, 10, 11, 12, 13, 14] have been of main interest 
in introducing the new hierarchy of the chaotic maps. On the other hand, the sensitively to 
the initial condition, control parameter and ergodicity which have tight relationships with the 
requirement of pseudo-random coding and cryptography [15, 16] are examples of interesting 
features of chaotic systems and it is natural idea to use chaos as a new source to construct 
new encryption systems [17]. 

In present paper, we introduce the rational order families of chaotic maps as a new hierarchy 
of chaotic map with an invariant measure. The Kolmogorov-Sinai entropy of these chaotic 
maps can be calculated analytically by using their invariant measure. An interesting property 
of these chaotic maps is their ability in simultaneous production and consumption of entropy. 
Additionally, being a measurable dynamically systems, so it can be studied analytically. 
The paper is organized as follows: In Section 2, we introduce the rational order families of 
chaotic maps. In Section 3, the invariant measure of these maps are given and in Section 4, 
we review the Kolmogorov-Sinai entropy and compute it for the rational order chaotic maps. 
Finally in Section 5 we calculate the Lyapunov exponent numerically and compare the results of 
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simulation with analytically calculated Kolmogorov-Sinai entropy The last sections contains 
our conclusion and two appendices. In these appendices we have calculated the invariant 
measure of the rational order families of chaotic maps via two different methods. 

2 Hierarchy of rational order families of chaotic maps 
with an invariant measure 

We first review hierarchy of one-parameter chaotic maps which can be used in the construction 
of families of rational order chaotic maps with an invariant measure. The one-parameter chaotic 
maps [1] are defined as the ratio of polynomials of degree N: 

a 2 (l + (-l)" 2 F ± (-N,N I*)) 



4>]f(x, a) 



(a 2 + l) + (a 2 -l)(-ir 2 F 1 (-N,N, 1 1 ,x) 
B?{T N {^)f 



2> 



- l + (a2-l)(T w (v^) 2 

(f) 2 N (x, a) 



a 2 (l-(-l)" 2Fl (-N,N,l(l-x))) 



( a 2 + 1} _ (a 2 _ 1)( _ 1} iv 2Fi( _^ N: i_ : (1 _ x)) 

B?{U N {^x))f 
1 + (a 2 - l)(U N (y/(l -x)f) 
where N is an integer greater than one. Also, 

2 F 1 (-N,N,^,x) = (-l) 7V cos(2iVarccos v ^) = (-1) N T 2N (^) 

is the hypergeometric polynomials of degree iV and T N (U n (x)) are Chebyshev polynomials of 
type I (type II), respectively. Here in this paper we are concerned with their conjugate maps 
which are defined as: 

4>n\x, a) = h o (/)pf(x, a) o h~ l = tan 2 (iVarctan y/x), 



4>${x, a) = h o 4$(x, a) o h 1 = 4 cot 2 (iV arctan -4=). 



(2-1) 

Conjugacy means that invertible map h(x) = ^f- maps / = [0, 1] into [0, oo). 

Now, in order to generalize the above hierarchy of integer order chaotic maps to the hierarchy 
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of rational order chaotic maps with an interesting property of simultaneous production and 
consumption of entropy, we need to replace x n+ i with a non-linear function of x n+ ±, particularly 
a non-linear function of chaotic maps of above type. But in order to have a single-valued map, 
we will take one of inverse branches of above nonlinear functions for x n+ i in each step with a 
probability equal to the probabilities of occurrence of the branches in iteration of the maps. As 
we will show in section 3, the probability of occurrence of each branch is equal to the integral 
of invariant measure of the map over the corresponding domain of the same branch. Therefore 
we can define these maps as: 

x n+1<k = g 2 ^ k o gi ( Xn ), «,j e{l,2}) with probability P k , (2-2) 

where functions g 1 and g 2 can be chosen as one the functions given in (2-1) and P k are proba- 
bilities of occurrence of inverse branches g^l of the map g 2 in its iteration. As we will see at 
the end of this section, the existence of an invariant measure will impose a relation between 
their parameters. 

Off course the functions given in (2-1) is not the only choice for the functions g\ and gi that 
leads to the hierarchy of rational order chaotic maps with an invariant measure. Obviously 
the following choices of the functions g\ and #2 



a) 


Tjr tan 2 (Narccoty/x), 


b) 


7p cot 2 (Narctany/x), 


c) 


7p cot 2 (Narccot^/x), 


d) 


^ tan(Narctan\x\)\, 


e) 


^ tan(Narccot \x ) , 


f) 


TjJ cot(Narctan\x\)\, 


g) 


;|| cot(Narccot\x\) , 


h) 


jt t an (Narctanx), 


I) 


7j tan( Narccotx), 


J) 


7j cot(Narctanx), 


k) 


7j cot (Narccotx) 







(2-3) 



lead to the hierarchy of rational order chaotic maps of trigonometric types (with an invariant 
measure), where some of them are equivalent to each others up to conjugacy. Also with the 
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choices of g\ and g 2 as [4] 



a) ^sc 2 (Nsc-\^)) b), ics^iVcs-^v^)) , 



(2-4) 



we get the Hierarchy of elliptic rational order chaotic maps of cs and sc types, where their 
invariant measure can be obtained for small enough values of module k of elliptic functions. 
Also it is possible to choose the function g\ and g 2 as one of the combined chaotic maps of 
Ref. [3]. 

Here in this paper we will consider the hierarchy of rational order maps with 

<7i(ai, Ni,x n ) = — tan(A r i arctana; n ) and g 2 (a 2 , N 2 , x n +i) = — tan(iV 2 arctan:r n+ i), (2-5) 
ai a2 

i.e., we have 



/ arctan(|^ tan(A 7 " 1 arctanx„)) k 2 n 
x n+lM = tan ^ ^-^ + — 



^ with probability Pk 2 , k 2 = 1, 2, N 2 , 

(2-6) 

where Ni and ^2 are integer greater than one and ai and a 2 are control parameters. As we are 
going to see in section 3, the maps (2-6) posses an invariant measure provided that we choose 
the parameters ai and a 2 in the form given in Equations (3-12) and (3-13), respectively. As 
an example we consider the following map for N x = 3 and N 2 = 2: 



ai 1 - 

x n+i,± — x -5 

a 2 3x n 



Q r 2 



ai 1 - 
1 + ( — x — 

a 2 3x r , 



3x 



n \ 2 



X v 



) 2 with probabilities P± 



(2-7) 



3 Invariant measure 



A dynamical system even time-discrete one-variable system has a number of possible types of 
behavior. The system can be in a fixed point and nothing changes, the trajectory of the system 
may also be on a cycle with a certain period. Fixed point and periodic orbits may be stable or 
unstable. We are usually interested in an invariant measure /x, i.e. a probability measure that 
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does not change under the dynamics. The probability measure \i on [0, 1] is an Sinai-Rulle- 
Bowen (SRB) measure as an invariant measure which describes statistically stationary states 
of system and absolutely continues with respect to Lebesgue measure. 

Now in order to determine the invariant measure of the analytical system described by the 
maps given in (2-6), we can write it as combination of the maps g\ and g^ 2 (as the fc 2 -the 
inverse branch of g 2 ) in the following form: 

x n +i,k 2 = 92M 9iM with probability P k2 k 2 = 1, N 2 (3-1) 

with g 1 and g 2 given in (2-5). 

Obviously the function g 2 (., a 2 , N 2 ) maps, its N 2 inverse branches x n+1 ^ 2 k 2 = 1,2, ...,N 2 
with corresponding different domains Ax n+ \ t k 2 (Ax n+ i^f) Ax n+ ij = for i 7^ j = 1,2, ...,N 2 ) 
into the same region. Therefore, if denote its value by y for different values of its argument 
then the map (3-1) can be written as: 

g 2 (x n+1 , a 2 , N 2 ) = y = g x (x n , a i? N^, (3-2) 

irrespective of to which branch or domain, the output x n+ i belongs (see Fig. 1). But in 
order to have a single output or single valued dynamical map, we have to consider only one 
of possible x n+ \ in each step with some probabilities or weights. Certainly the most natural 
weight of a given branch is the corresponding probability of its occurrence in infinite iteration 
of map y — g 2 (.,a 2 , N 2 ), where it can written in terms of its invariant measure \i 92 as, 

P(occurrence of k 2 — the branch) = / /j, g2 (x)dx. (3-3) 

Therefore the invariant measure of this map should satisfy the following Frobenius-Perron 
integral equations: 

fi(y)=f 5(y - g 1 (x n ,a 1 ,N 1 ))ii(x n )dx n (3-4) 
J 

and 

Piv) = I ~ 92(x n+1 ,a 2 ))fi(x n+1 )dx n+1 , (3-5) 
J 
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which are equivalent to: 

v(y) = 

and 



where 



and 



v(y)= 2 M g n+i,fc 2 )l <fa: ^ / 1 ' fc2 l» ( 3 - 7 ) 

^n+l,fc 2 G9^ (2/) 

( /c^ 7T 

x n>fcl = tan(— arctan(aiy) + — ), fc 2 = 1, N v 
x n+1M = tan(-^- arctan(a 2 y) + ^), k 2 = 1, iV 2 . 

iV 2 iV 2 

The invariant measure /i(y) for <7i(x) can be written as: 

/ x ^ a, / l + x 2 nkl \ , , 

Assuming that has the following form: 

= WVW) < 3 " 9 » 

where for (5 — 1, it reduces to the invariant measure which has already applied to pushout 
measure [18], expression (3-6) reduces to 

\2 Nt 



By comparing of both sides of Equation (3-8), we can determine a^, (i = 1, 2) as 

a, = , (3-11) 

2^fc=0 °2fc+lP 

for even values of N i: and 

a, = ° 2fc +^ , (3-12) 

for odd values of iVj,(for proof see Appendix A). 

Therefore ai and a 2 depend on the parameter (3 and integers N\ and iV 2 , respectively. Also to 
make the paper more readable, we have derived the invariant measure of the map 
y = \ tan(4arc£cmx) by using Shure's invariant polynomials in Appendix B. 
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4 Kolmogorov- Sinai entropy 

In this section we review first, the Shannon entropy and then talk about Kolmogorov-Sinai 
entropy (for more details see [19]). Consider dynamical system characterized by a certain 
iterative map. Let B = (Bi, Bj, . . . , B n ) be a decomposition of the unit interval along x n . 
Now we subdivide each interval Bi into say A points, and perform ( iterations on each one 
of them so we make sure that transients have died out. Then A points by then will spread 
to other subintervals. A percentage of them will be perhaps located within the limits of Bj . 
After transients die out the common area of F^(Bi) and Bj, e.g. F^(Bi) f]Bj will be express 
in a non-normalized way the number of elements of Bi reaching Bj after £ iterations. So in 
normalized form: 

wHBjlBi) = ^y> , m 

here //(.) = J c /i(x)dx, where c is the pertinent interval. The entropy of the chosen partition 
or, the average amount of information needed to locate the system in state space is given by 
the Shannonian entropy; 

A 

S=-E ^ log 2 vWbits. (4-2) 

i 

The A values fi(Bi) may be calculated from the Wij elements from the (A — 1) equations of 
the linear system: 

M^ = EM« (4-3) 

and the normalization condition: 

!>(£,■) = i, 

where the transition probability matrix describes the probability of jumping in one step 
(iteration) from the element Bi of the partition to the element Bj. The average amount of 
information created by the linguistic system by per transition per unit time is given by the 
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Kolmogorov-Sinai entropy for the chosen partition; namely; 

A A 

Sk = E E vWWij log 2 Habits. 

i=l j=l 



(4-4) 



The macroparameter however, characterizing the degree of grammatical coherence of the cre- 
ated Markovian chain is the mutual information or transinformation. 



m = ee ti^m n b s ) i og2 f,l ^ in . i)nn /' \^. 

i=i i=i 



(4-5) 



It stands for the information stored in a symbol along the sequence about what is going to 
emerge £ iterations (or £ time units ) later, /(£) gives the information transferred between two 
symbol £ steps apart. As the number of the decomposition of the unit interval goes to infinity, 
in such a way that, the size of each intervals {Bj) goes to zero. The mutual entropy given in 
(4-5) reduces to the well known Kolmogorov-Sinai (KS) entropy which is given by: 



d 

fi(x)dx In | —g(x, a, N) \ 



(4-6) 



with g(x, a, N) = tt (tan( iVarctanx)) h(fx, g(x, a, N)) can be written as: 



/+oo 
-77 
-oo 7T 1 



VP 



(4-7) 



-oo n(l + (3x 2 ) 1 a 1 +a; 2 

Following the calculating of Ref.fl], one can show that after a change of variable Vfix = tan 9, 
and using the integral of type; 

j n I a+\/a 2 -b 2 I 



1 / >7r 

— / In a + 6cos#| = < 
7T Jo 



In 



|a| > |6|, 
l a l — 1^1) 



(4-8) 



we get the following expression the KS-entropy: 



h(fi,g(x,&,N)) = -j=ln 



N { v / ^(EU'o 1 ^^)a + a(EL 2 o 

Q 3 



(4-9) 



(v^+l)(ES^V) 

Now, we come to calculate the KS-entropy of fractional order maps. Before getting to involved 
with the details of calculation, we first talk about simultaneous production and consumption 
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of entropy in these maps. Figure 1, gives us an insight to see how this is possible. In this 
figure, Ni ( N 2 ) corresponds to the number of branches x n+ i( x n ) of the map at time n+ l(n). 
The left hand half of figure 1 shows the contraction of Ni identical branches of x n+i into 
a single branch, while the right hand half of the figure shows branch out to N 2 branches 
x n , where they correspond to increment and decrease of entropy, respectively. It should be 
reminded that in contraction of branches, entropy increase due to loss of information, while in 
branch out it decreases due to reception of information. Now in order to calculate KS-entropy 
of fractional order maps(2-6), we should notice that the right and left halves of Figure 1, 
correspond to the right-hand and left-hand sides of Equation (3-3) ( equivalent of (2-6)). In 
other words, there are Ni convergent branch x n+ i and N 2 divergent branch x n in the left and 
right halves of Figure 1. Therefore, according to Figure 1 there are N 2 possible final states 
x n+ i t k 2 , k 2 = 1,2, ...,N — 2 with corresponding weights given in (3-3), where we should take 
average over their corresponding KS-entropy given in (4-6). Hence, the KS-entropy of rational 
map (2-6) can be written as 



r r+oo 

h(/i, rational order map) = / /i(x n+ i, a 2 )dx n+ i / n{x n , &i)dx n In ' 

fe 2 = l •* ^ x n+l,k 2 J-OO 

= fi(x n+1 , a, 2 )dx n+1 / u.(x n ,a,i)dx n W 

U 1 JAln+li, J-OO 



n+l 



= Yj! n{x n+1 ,aa)dx n +\ f fi(x n , ai)dx„(ln | ^ n+1 1 + In |— —\) 

h i •> &x n +i,k 2 j-oo ay ax n 

/+°° dy f +oc dy 
//(x n ,ai)da; n ln|- — \- fi(x n+1 , a 2 )dx n+1 \n\- 1, (4-10) 
-oo ClX n J — oo dX n -\-i 

where in the last line above we have used the following normalization relation 

/+0O 
u,(x, cLi)dx =1 i — 1,2. 
-oo 

Now, comparing the last line of (4-10) with (4-6), we get the following expression for KS- 
entropy of fractional order map (2-6): 

/i(rational order map) = h(/i, — tan(Aq arctanx„)) — h(/i, — tan(A r 2 arctanx n ) (4-11) 

Hi a 2 
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with h(/i, gj-tan(iVjarctan£„)) i — 1,2 given in(4-9). 

Obviously formula (4-11) implies the simultaneous production and consumption of the en- 
tropy, where the term with positive sign corresponds to the production of the entropy while 
the term with minus sign corresponds to the consumption of the entropy, respectively. Also, 
it is interesting to note that maximum value of the entropy is equal to ln 2 that corresponds 
to, ai = a 2 = 1 ( actually this is the main reason for naming these maps as rational order maps). 



5 Lyapunov exponent and simulation: 

A useful numerical way to characterize chaotic phenomena in dynamic systems is by means of 
the Lyapunov exponents that describe the separation rate of systems whose initial conditions 
differ by a small perturbation. Suppose that there is a small change 5x(0) in the initial state 
x(0). At step or time n this has changed to Sx(n) given by: 

5x(n) » MO) 1^1 = 5x(0)\^^.p^ (5-1) 

dx dx n -i dx n -2 dxo 

where we have used the chain rule to expand the derivative of ^ 2L . In the limit of infinites- 

r dxo 

imal perturbations 8x(0) and infinite time we get an average exponential amplification, the 
Lyapunov exponent A, 

A = limn^oo - In ^ X ™ | = lirrin^-ln |^^| = lim n ^<J- V In | f Xk |. (5-2) 
n dx{0) n dxo n ^ dxu-i 

Similarly the Lyapunov exponent of rational maps (2-6) can be obtained from formula (5-2) 

provided that we replace Xk, at random with Xk,± at each step k with probabilities P± — \. 

Here in this work we have simulated Lyapunov exponent of rational order map with A r 1 = 2 

and N 2 = 3 for different values of j3, where the result that obtained supports the analytic 

calculation of KS-entropy given in (4-11) for particular case of N± — 2 and N 2 = 3 (see Figure 

2). 
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We have given a new hierarchy of rational order families of chaotic maps which presents an 
interesting description of simultaneous production and consumption of entropy. Using the SRB 
measure, the Kolmogorov- Sinai entropy of chaotic maps have been calculated. It would be 
interesting to introduce these kinds of maps in higher dimensions which is under investigation. 



7 Appendix A 

Here in this appendix following the prescription of References [1, 2, 3, 4, 5], we prove that 
the invariant measure given in (3-9) satisfies the corresponding PF equations of the maps 
Hi = gj- tan(iVj arctanxj), i — 1,2, provided that the parameters ai and a 2 can be expressed 
in term of (3 as in formulas (3-11) and (3-12). To do so we we can write the right hand side of 
Equation (3-8) as, 

where we have omitted the indices i and n. Hence, Equation (3-8) can be written as: 

1 + a V _ a a(/?-l) d ln( yr (3 -i 2)) (? 2) 

To evaluate the second term in the right hand side of above formulas we can write the equation 
V = \ tan(iV arctanx) in the following form: 

= aycos(A^arctanx) — sin(iV arctanx) 



f] [VI \ 



—^EUvY,Cl{-l) k x 2k -x E C 2 N k+1 (-l) k x 
constant . . 

(1+X 2 )2 k=1 
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where Xk = tan(-^ arctan(ay) + ^) k — 1, N are its roots. Therefore, we have: 

^f=i in (n (r 1 + 4)) = (n + **)) + in (n + **)) ) 



9/5- 1 



In 



9/3- 



ln 



(1 - /T 1 )^ ^aycos(iVarctan(-iy / /Fi)) - sin(N arctan(-iy / /F*)) 
(1 - /T 1 )^ ('aycos(Narctan(iy / /Fi)) - sin(ATarctan(iy / /Fi)) 
9 



dp~ 



ln 



ay(E01 2 +r'(E01 2 



fc=0 



fc=0 



d 



ln 



(1 - P~Y [afy 1 cos\N axctsm(i^ (3- 1 )) - sin 2 (iV arctan^/?- 1 )) 
7V/3 /3iV(l + aV)A(l)i?(i) 



+ 



z^- 1 (i-nMS))w+(B(j)) 2 ) 

with polynomials A(x) and defined as: 

[f] 



(7-3) 



fc=0 



In deriving the of above formula we have used the following identities: 



cos(iV arctan i\fx) 



A(x) 



sin (N arctan iyfx) = i\fx 



(1-x) 
B(x 



N_ 1 

2 



(1-x) 



N_ 1 

2 



inserting the results (7 — 3) in (7 — 2), we get: 

1 + aV 1 + aV 



i + Py 2 



g(j) . ai &A ^l\ 2 



Hence to get the final result we have to choose the parameter a as: 



a = 



4*) ' 



(7-4) 



(7-5) 
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8 Appendix B 

Here we derive the invariant measure of chaotic maps by using Shure's invariant polynomials 
in a way which is different form that of Appendix A. In order to make the paper more readable 
we consider only N — 4 case, i.e, the map,. 

y = -tan(4arctan:r) (8-1) 
a 



which can be written as: 



or 



4x(l-x 2 ) . , 

ay = T^rh (8 - 2) 



ay ay 

The needed Shure's invariant polynomials of variables xi,x 2 , ...,X4 are defined as: 

Si = x 1 + x 2 + x 3 + x 4 , 

S U = XiX 2 + X±X 3 + X\X 4 + X2X3 + X 2 X 4 + X3X4, 

Sin = x x x 2 x 3 + x x x 2 X4 + xxx 3 x 4 + x 2 x 3 x 4 , 

Sim = xxx 2 x 3 X4 (8-4) 

and 

rr 2 , 2 , 2 , 2 

o 2 — X-y "T X 2 "T Xq ~r X 4 , 

q 22, 22i 22i 22i 2 2, 22 

O22 — ^1*^2 I *^1^3 I X -^X ^ I 3'2^'3 I Xc^X ^ I 3/2^4.5 

O 222, 222, 222, 222 

w <-y»'-' _l_ .-y* /"v* ■"Y* —I— T* T* r f* —I— 'T* <"y» <"y» 

^222 — 1 2 3 1 2 4 13 4 ^ ^' 2 3 4 ? 

Sim = x\x\x\x\. (8-5) 
The second set of Shur's invariant polynomials can be expressed in terms of the first one as: 



Rational order chaotic maps 16 

S22 = 5*11 — 2S'iS'iii + 2S1111, 
5*222 — 5*111 ~~ 2<S'ii<S'iiii, 

^2222 = S'lin- (8-6) 

Considering the variables xi,x 2 ,xs and rr 4 as roots of Equation (8-3), one can obtain the first 

set of Shur's invariant polynomials given in (8-4) as: 

4 4 

Si = , Sn = —6, Sui = — , Sim = 1 (8-7) 

ay ay 

and using Equation (8-6), we have: 

16 32 16 

5*2 = 7 — + 12, S 22 = 7 — vj + 38, 5222 = 7 — + 12, S 2222 = 1- (8-8) 
(ay) 2 (ay) 2 (ay) 2 

Again, writing the PF equation of map (8-1) and assuming that its invariant measure is of the 

form (8-3), we have: 

The summation on the right-hand side can be written as : 

« 1 + xj (1 + x\){\ + (3xp(l + (3xp(l + (3xt) rRim 

ki + P4 nti(i + /^) ' 1 j 

where using Equation (8-8), the numerator of above fraction becomes: 

4 + (3/3 + 1)S 2 + (2/3 2 + 2[3)S 22 + (/3 3 + 3f3 2 )S 222 + 4[3 3 S 2222 = 16(1 + (3){l + 6/3 + /3 2 ). (8-11) 
Also using again (8-8), for its denominator 

1 + (3S 2 + /3 2 S 2 2 + P 3 S 222 + [3 A S 2222 , (8-12) 

we get: 

16 /3(l + /3 2 ) + (/3 2 + 6/3 + l) 2 . (8-13) 



a 2 y 2 



Using the results that obtained above, the invariant measure takes the following form : 

„( V ) = Ml + ^)(l + 6/3 + ^ 2 ) 

^ aV(/3 2 + 6/3 + l) 2 + 16/3(/3 + l) 2 1 ; 
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(8-15) 



which should be equal to 1+ ^ 2 , where this is possible only for the following choice of a 

4(3 + 4(3 2 
a ~ 1 + 6(3 + (3 2 
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Figures captions 

Fig. 1. The schematic diagram of forward and backward branching of rational order chaotic 
map for describing simultaneous production and consumption of entropy. 
Fig. 2. Shows Lyapunov exponent (solid curve) and KS entropy (o) versus the control param- 
eter (3, as this figure shows, there is a maximum at (5 which corresponds to a x = a 2 = 1. 
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